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The unique forward LHCb acceptance opens the possibility for performing precision Standard 
Model (SM) processes in a kinematical regime beyond the reach of both ATLAS and CMS exper¬ 
iments. In this article we discuss the feasibility of performing W ± j measurements with 14 TeV 
LHCb data, and find that the leptonic asymmetry at high pseudorapidity, with appropriate cuts 
on the associated jet, can provide a factor of three reduction on the d quark PDF uncertainty in a 
region of Bjorken x > 0.5. 


INTRODUCTION 

With the large amount of data expected in Run II-A 
and Run II-B at the LHC, corresponding to pp collisions 
at a centre-of-mass (COM) energy of 13 and 14 TeV re¬ 
spectively, precision Standard Model (SM) measurements 
and beyond-SM searches will become possible in kine¬ 
matic regimes which have so far been statistically inac¬ 
cessible. In both cases, it is necessary to reduce theoret¬ 
ical uncertainties to improve the sensitivity of searches 
and precision SM measurements. 

One of the major uncertainties for many measurements 
at the LHC are those associated to the modelling of ini¬ 
tial state partons. Fully describing the inelastic scatter¬ 
ing of two hadrons is complicated due to the fact that 
the colliding hadrons are strongly-bound composite ob¬ 
jects of interacting quarks and gluons. This complication 
has previously been overcome by factorising the full scat¬ 
tering process into both short-distance and long-distance 
effects PQ . The short-distance effect describes a partonic 
interaction occurring at large momentum exchange scale 
Q while the long-distance effects describe the evolution 
of partons within the incoming protons from a low en¬ 
ergy scale Q 0 to the scale of the short-distance interac¬ 
tion Q. Although the partonic evolution can be com¬ 
puted perturbatively, as Q o is chosen to be well above 
the scale Aqcd at which the strong interactions are non- 
perturbative, the distribution of the partons within the 
hadron at Q o is a priori unknown. Rather, these par- 
ton distribution functions (PDFs) are extracted exper¬ 
imentally, based on a global analysis of data of many 
different measurements performed at many different en¬ 
ergy scales. Several groups now provide modern PDF 
sets obtained from global fits to a large sample of data 
sets: ABM (2]; CTEQ-TEA [3]; CTEQ-JLAB g]; HER- 
APDF [5]; MMHT |J| andNNPDF [7j. However, beyond- 
SM searches and precision measurements of SM processes 
are often limited by PDF uncertainties associated to the 
extraction of the PDFs from the global analysis. Im¬ 


proving the constraints on these PDFs will enhance our 
ability to test the SM and to look for new physics signals 
which may be kinematically accessible in Run II. 

Due to the forward acceptance of LHCb, measurements 
of SM processes naturally provide information on the 
colliding partons which simultaneously carry a relatively 
large and small momentum fraction of the incoming pro¬ 
tons momentum. As the probability of finding a parton 
declines with increasing x, measurements in which both 
partons are at large-a; are less likely than those with only 
one parton at large-a;. Consequently, measurements at 
LHCb provide the best environment for constraining the 
large-a; region of PDFs. Indeed, several recent fits already 
include Run I LHCb data from various electroweak mea¬ 
surements Ml, which has improved the uncertainties of 
the u and d distributions at large-ax In this paper we sug¬ 
gest that measurements of W ± production in association 
with jets within the LHCb acceptance can be used to con¬ 
strain these PDFs at extremely large-a;. In this letter we 
discuss the feasibility of W ± j cross section and asymme¬ 
try measurements with Run II-B data and the expected 
improvement in the d-quark PDF at large-a;. Further to 
this, we also demonstrate that requiring a significant pr 
cut on the reconstructed jet provides sensitivity to values 
x > 0.5. Such data is of general importance to global fits 
as the available data in this region from deep-inelastic 
lepton-hadron scattering is statistically limited. 


CROSS SECTION PREDICTIONS 

The goal of this article is to evaluate the potential sen¬ 
sitivity of Wj measurements at LHCb to PDFs at large- 
ax The presence of a reconstructed jet in the final state 
acts to increase the partonic COM energy, predominantly 
providing sensitivity to the PDFs of proton travelling in 
the direction of the LHCb detector at large-ax A signifi¬ 
cant pp cut on the reconstructed jet is therefore vital for 
improving the sensitivity to PDFs. To avoid otherwise 
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overwhelming QCD backgrounds, only the leptonic final 
state pp —> Ivi j is considered, and the provided differ¬ 
ential cross sections correspond to only a single lepton 
flavour final state. Importantly, as the LHCb acceptance 
is limited to the forward-region, fully reconstructing the 
W boson using transverse momentum conservation is not 
possible, due the ambiguity that missing transverse mo¬ 
mentum may be due to undetected neutrinos or visible 
particles which are not within the detector acceptance. 

The generation of the IVj events is performed with 
MadGraph5_aMC@NL0 m subsequently showered with 
Pythia8 m , achieving NLO+LL accuracy. The NLO 
NNPDF3.0 PDFs a s (m z ) = 0.118 0 are used for the 
generation of the hard process, and the default Pythia8 
Monash. 2013 tune m is used. As constraints at large-a; 
are also of interest to the fitting of nuclear corrections to 
lepton- or proton-deuteron scattering processes, we also 
provide comparisons using the NLO CJ12 PDFs [J]. 

Charged leptons are required to have a pr greater 
than 20 GeV, and to be within the LHCb acceptance 
2.0 < rji < 4.5. To coincide with previous LHCb mea¬ 
surements which involve jet reconstruction m , jet recon¬ 
struction is performed using the FASTJET3 software [T5] 
with the anti-fct algorithm El and a distance parameter 
of R = 0.5. The jets are also required to be within the ac¬ 
ceptance 2.0 < 77j et < 4.5. In the event selection, at least 
one reconstructed jet within the acceptance must have a 
Pt greater than p™* (which we vary), and an isolation cri¬ 
terion is placed on the lepton such that Al?(/ ± Jet) > R. 
Sub-leading jets are only reconstructed if they have a pr 
greater than 10 GeV. 

In Figure [l] the differential cross section with respect 
to the lepton pseudorapidity is shown for pf/ ut = 60 and 
80 GeV. A scale uncertainty is computed by varying fac¬ 
torisation (pf) and renormalisation (pr) sales indepen¬ 
dently by a factor of two around the nominal scale ( po ) 
with the constraint 1/2 < Pf/pr < 2. The nominal 
scale is set dynamically on an event-by-event basis to 
po = (XZ, Pt)/2, where i corresponds to the W boson 
and jets. The shown PDF uncertainty is the 68% confi¬ 
dence interval. It is worth noting that the jet passing the 
p™* within the LHCb acceptance is the leading jet for ~ 
0.95 events, suggesting that the generation of W ± +lj at 
NLO is sufficient for the considered final state. The scale 
uncertainties (6-16)% are dominant across the considered 
phase space, however the PDF uncertainties become sig¬ 
nificant for large values of lepton pseudorapidity. For 
the extreme lepton pseudorapidity bins (3.9 < pi < 4.5), 
the average < X\ > values are 0.40, and 0.43 for p™* = 
60, 80 GeV respectively. For the considered PDF set, the 
shown uncertainty predominantly arises from lack of con¬ 
straints on the d-quark PDF for large-a; and Q 2 values. 

In Figure [2] the normalised differential cross section 
for the W~j subprocess is shown. The distribution is 
normalised with respect to the differential distribution 
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FIG. 1. Differential cross section with respect to charged 
lepton pseudorapidity. The sum of W^j events are considered, 
and a single leptonic branching fraction has been applied. 
Distributions are shown for pfjV values of 60 (blue) and 80 
(red) GeV where scale (solid) and PDF (dashed) uncertainties 
are displayed. For clarity, the relative uncertainties for the 
two sources are shown separately below the figure. 

integrated over the experimentally accessible region cor¬ 
responding to pi £ [2.0 —4.5] — referred to as the fiducial 
cross section. As less stringent constraints are present for 
the d-quark distribution compared to the u-quark, we 
have chosen to present only the W~j subprocess as the 
corresponding PDF uncertainties are largest. By normal¬ 
ising the distribution, the renormalisation scale uncer¬ 
tainty is reduced as the leading contribution to this un¬ 
certainty, through the running of a s (Q 2 ), partially can¬ 
cels for this observable. For the pseudorapidity range 
p £ [2.6 — 3.3], the average Q 2 of these bins is similar 
to the fiducial cross section which normalises the distri¬ 
bution, and causes the partial cancellation. The PDF 
uncertainties are less sensitive to this normalisation in 
the high and low-pseudorapidity bins as they depend on 
both x and Q 2 , and the x dependence of the fiducial cross 
section and the high/low pseudorapidity bins is less cor¬ 
related than the Q 2 dependence. 

ASYMMETRY PREDICTIONS 

Beyond normalised differential cross sections, a mea¬ 
surement of the charge asymmetry as a function of the 
lepton pseudorapidity is extremely well motivated. This 
observable is robust with respect to higher-order correc¬ 
tions as they do not differentiate between W + j and W~j 
subprocesses. The scale uncertainties only contribute 
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FIG. 2. Normalised differential cross section for the W ~j 
subprocess with respect to charged lepton pseudorapidity for 
Px ut values of 60 (blue) and 80 (red) GeV. The cross section 
is normalised with respect to the fiducial cross section. Both 
scale (solid) and PDF (dashed) uncertainties are displayed. 


through factorisation dependence, which can affect the 
w/d-quark ratio or the quark-antiquark PDF at a given 
value of x and Q 2 , both of which alter the asymmetry. 
Experimentally, this observable is also well motivated 
since many systematic uncertainties cancel. While not 
explicitly considered as a source of uncertainty, it has 
been shown 033 that the choice of parton shower tune 
can also have an effect on the jet px distribution. Asym¬ 
metry predictions obtained alternatively using the older 
Tunel and Tune4Cx tunes of Pythia8 show deviations 
of the asymmetry predictions on the order of the scale 
uncertainties. 

In Figure [3j the differential asymmetry is shown for 
the choice p™ 1 = 80 GeV. The expected statistical pre¬ 
cision after 5 fb _1 of data taking is also shown and is 
obtained by scaling the expected yield of events by the re¬ 
construction efficiencies reported in existing LHCb mea¬ 
surements. In particular, muon and jet reconstruction 
efficiencies of 0.75 and 0.95 respectively are determined 
in measurements of inclusive W ± production [18 and 
Z production in association with jets [14] and are used 
in this study. In addition to the NLO+PS prediction, 
the corresponding central fixed-order prediction is also 
included showing satisfactory agreement. 

Noticeably, the lepton asymmetry tends to decrease 
beyond rji > 2.6. Although the W boson asymmetry 
increases at large rapidities, as the zt/d-quark ratio and 
the quark-antiquark PDF asymmetry increase at large- 
x, the fraction of W bosons produced in a left-handed 
polarisation also increases mm- As a consequence, the 
lepton asymmetry is reduced as the charged lepton from 
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FIG. 3. The charged lepton asymmetry for selected events 
containing a jet with a pt of greater than 80 GeV. Both scale 
(yellow) and PDF (green) uncertainties are shown, as well 
as the projected statistical uncertainty on measurements at 
LHCb with 5 fb _1 of data. The mean value of x\ is shown 
below with the band representing the root mean square of the 
distribution. 

the polarised W + ' K ~' > boson decay prefers to be produced 
at lower(higher) rapidity with respect to the W boson 
rapidity. 

The relative difference of the predicted asymmetry 
with respect to NNPDF is shown for both CJ12min and 
CJ12max PDF sets in Figure [4] The statistical uncer¬ 
tainties from Figure [3] are now centred on the CJ12max 
prediction. For CJ12max, the PDF uncertainty is explic¬ 
itly shown and corresponds to the ler Hessian uncertainty 
with a tolerance factor, T = 10. While not explicitly dis¬ 
played, a similar size uncertainty is also found for the 
CJ12min set. At high pseudorapidity the CJ12min and 
CJ12max predictions become differentiable. The differ¬ 
ence between the shown CJ12 predictions is the assump¬ 
tion on the strength of nuclear corrections employed in 
the simultaneous fit of PDFs and nuclear corrections (re¬ 
quired to model the deuteron bound state) in the CJ12 
global analysis, which suggests that LHCb data can indi¬ 
rectly provide useful information on the strength of these 
corrections when included in such a fit. 


ESTIMATED PDF UNCERTAINTY REDUCTION 

To quantify the sensitivity of the proposed asymme¬ 
try measurement, a reweighting of the NNPDF3.0 NLO 
replica PDFs is performed following the procedure out¬ 
lined in Emu- Although no LHCb data is currently 
available for this observable, the reweighting procedure 
can be performed with respect to the central PDF predic- 
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FIG. 4. The relative difference of the predicted charged lepton 
asymmetry with respect to the NNPDF prediction of Figure 
[3] is shown for the CJ12min (red) and CJ12max(blue) PDF 
set. The PDF uncertainty is shown for the CJ12max PDF set. 
The statistical precision of the LHCb dataset with 5 ffcff 1 is 
also shown, where the central value is fixed to that of the 
CJ12max prediction. 
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tion. This allows the reduction of the relative PDF un¬ 
certainty to be quantified when assuming different levels 
of experimental precision. 

This precision is defined as the relative uncertainty 
on the ratio of W + j to W~j production, R W j , which 
is related to the asymmetry, A, through the formula 
A = (Rwj — 1 )/(Rwj + !)• While the relative uncer¬ 
tainty on A depends on the central value and is largest 
for asymmetries close to zero, the relative uncertainty on 
R Wj provides a more meaningful estimate of the exper¬ 
imental precision, and can be translated into an uncer¬ 
tainty on the asymmetry for any given central value. 

In the measurement of inclusive W production at 
LHCb [18] , the uncertainty on the ratio is dominated 
by the purity determination and is reported to vary be¬ 
tween 1 and 3% across the different pseudorapidity bins. 
The addition of a jet to the final state is not expected 
to inflate the uncertainty to any great extent, as the un¬ 
certainties related to jet reconstruction will also largely 
cancel in the ratio. While the finite jet energy resolution 
of the detector, reported to vary between 10-15% in jT4] . 
will result in a migrations of events above and below the 
chosen jet px threshold, a smearing of the jet energy by 
this amount in the simulated samples indicates that the 
effect will be small. Consequently, the precision of the 
measurement is expected to be at the level of a few per¬ 
cent. 

The reweighting is performed assuming 5/50 fb -1 of 
data is collected, and the impact on PDFs is studied if 
the additional systematic uncertainties can be controlled 
at the level of 1 or 2%, which are added in quadrature 
with the statistical uncertainty in each bin. The sample 
with a pt threshold for the leading jet of 80 GeV (see 


Fig.§, which contains eight independent data points, is 
considered. To begin with, the probability distribution 
for the PDFs is provided by a finite ensemble of 100 repli¬ 
cas. The potential impact of this data is then estimated 
by computing a weight (wj) for each of these replicas as 
detailed in m, which is conditional on both the old data 
as well as the new data points which we are including. 
After reweighting the number of effective replicas can be 
quantified by using the Shannon entropy as: 


N eS , = exp 


— Wj Log (N/vjj ) 


(1) 


After reweighting assuming 50 fb -1 of data, the original 
100 replicas are reduced to 18 and 30 effective replicas 
when a systematic uncertainty of 1 and 2% are assumed 
respectively. When performing reweighting with respect 
to the central value (as is done here), a reduction in the 
number of effective replicas indicates that the new data is 
highly constraining. The relative 68% confidence interval 
of the original and reweighted evolved d-quark PDF are 
shown in Fig. [5] 
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FIG. 5. Relative lcr uncertainty of the evolved d-quark PDF. 
The reweighting is performed with respect to the central PDF 
prediction, and the sensitivity of different experimental pre¬ 
cisions are also included. 


The d-quark PDF is already well constrained below 
x = 0.5, however, significant improvements can be ex¬ 
pected beyond this with the addition of LHCb data. As 
indicated in the lower plot of Fig. |3j although the RMS 
of sampled x values for high pseudorapidity events lies 
within the region x £ [0.35,0.60], a significant fraction 
of events contribute both above and below this region. 
The observed PDF constraints for NNPDF3.0 are most 
significant at large-x values. For example, at x = 0.7, 
assuming a 1% systematic and 50 fb -1 of data, the rela¬ 
tive reduction of the uncertainty of the reweighted with 
respect to unweighted PDF is ~ 60%, corresponding to 
a factor of three reduction. Even when 5 fb -1 of data 
is collected (expected after 1 year of data taking), the 
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corresponding reduction is ~ 35% indicating that early 
Run-II measurements will also be highly constraining. 

It is worth noting that although the improvement in 
the PDFs has been quantified by only considering the 
total uncertainty on the d-quark PDF at Q 2 = M 2 V ; 
correlated improvement in the other PDFs, which has 
not been quantified here, is also expected. For example, 
a reduction in the uncertainty for the u-quark PDF and 
d/u -quark PDF ratio at large-a; is expected, as well as for 
the PDFs at low-a; which also contribute the hard pro¬ 
cess. The substantial reduction in the number of effective 
replicas in this analysis also suggests that a new global 
fit will be necessary to fully quantify the impact of this 
data when it becomes available. This can be seen explic¬ 
itly in Fig. [ 5 ] in the region x £ [0.4, 0.5] where a marginal 
increase in the d-quark PDF is observed for Q 2 = m\ v . 
This demonstrates that the weighted ensemble provides 
a poorer description of the underlying probability distri¬ 
bution. 


DISCUSSION AND CONCLUSIONS 

We have proposed that measurements of the Wj pro¬ 
cess with a pr threshold for the leading jet of ~ 60, 
80 GeV should be performed at LHCb. Of particular 
interest is the lepton asymmetry which is both theoret¬ 
ically and experimentally clean, and sensitive to PDFs 
beyond x = 0.5. In this work, we focussed on the poten¬ 
tial of the Wj final state without taking particular care 
to optimise kinematic cuts. We note that placing further 
restrictions on the final state, such as requiring the lead¬ 
ing jet to also be forward (rjj > 3.0), or requiring the 
presence of a secondary jet within the acceptance with a 
moderate pt threshold (p™* > 40 GeV) further improves 
the sensitivity to yet larger values of x. Given that the 
proposed measurement is not anticipated to be statisti¬ 
cally limited, these additional requirements should also 
be considered in future measurements. For example, a 
measurement of the lepton asymmetry differentially in r/i 
as well as jet pp and/or rjj. 

The proposed measurements have the potential to pro¬ 
vide among the most stringent large-a; constraints from 
Run-II of the LHC. This is a consequence of the theo¬ 
retical and experimental precision of the chosen variable, 
as well as the unique forward acceptance of the LHCb 
detector, which naturally provides sensitivity to large-® 
values as compared to ATLAS or CMS detectors. Recent 
NNLO calculations .23' will also increase the sensitivity 
of LHCb measurements to the large-® gluon, which can 
then be used as a constraint as outlined in Ref. [23! . 

Precise differential Wj cross section and asymmetry 
measurements are also important as a background to 
other interesting SM measurements which can be per¬ 
formed in the LHCb acceptance. For example, in other 
work [2oIB71 . it has been suggested that top-quark pair 


cross section and angular production asymmetry mea¬ 
surements are possible at LHCb. In both cases, the 
statistically most promising final state is tt —> IbX , for 
which the Wj subprocess is the dominant background, 
and consequently measurements of Wj production can 
significantly improve the precision of tt measurements at 
LHCb. 
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